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Transfer-matrix methods are used to calculate spin-spin correlation functions (G), Helmholtz free 
energies (/) and magnetizations (m) in the two-dimensional random-field Ising model close to the 
zero- field bulk critical temperature T c q, on long strips of width L — 3 — 18 sites, for binary field 
distributions. Analysis of the probability distributions of G for varying spin-spin distances R shows 
that describing the decay of their averaged values by effective correlation lengths is a valid procedure 
only for not very large 7?. Connections between field- and correlation function distributions at high 
temperatures are established, yielding approximate analytical expressions for the latter, which are 
used for computation of the corresponding structure factor. It is shown that, for fixed R/L, the 
fractional widths of correlation-function distributions saturate asymptotically with L~ 2 2 . Consider- 
ing an added uniform applied field h, a connection between f(h), m(h), the Gibbs free energy g(m) 
and the distribution function for the uniform magnetization in zero uniform field, Po(m), is derived 
and first illustrated for pure systems, and then applied for non-zero random field. From finite-size 
scaling and crossover arguments, coupled with numerical data, it is found that the width of Po(m) 
varies against (non- vanishing, but small) random-field intensity Ho as Hg 3 ^ 7 . 
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I. INTRODUCTION 

The random-field Ising model (RFIM) has posed a 
number of challenges to researchers since its introduction 
as an apparently purely theoretical puzzle jlj . The later 
realization that it corresponds, give or take a few (hope- 
fully irrelevant) details, to the experimentally realizable 
dilute Ising antiferromagnet in a uniform applied field (^] 
brought new insights, and new questions as well; among 
the latter, is the interpretation of experimental data in 
a suitable theoretical framework. This has proved to be 
rather an intricate subject, even down to basic aspects 
such as whether the lower critical dimensionality for the 
problem was d — 2 or 3 || || . Though by now this partic- 
ular issue has been settled in favour of d — 2 Q , several 
important aspects (such as the scaling behaviour near 
the destroyed phase transition in d = 2, which will be of 
interest here) still require further elucidation j7| . 

In the present paper we deal with the two-dimensional 
RFIM, where long-range order is destroyed, and a zero- 
temperature, zero-field "anomalous" critical point ap- 
pears ||. The latter will not concern us directly, as we 
shall be working at high temperatures, close to the pure- 
system ferro-paramagnetic transition. We extend and 
complement our early work Q, making use of transfer- 
matrix (TM) methods on long, finite-width strips of a 
square lattice; wc generate and analyze statistics of spin- 
spin correlation functions and uniform magnetizations. 
Wherever feasible, we attempt to draw connections be- 
tween our numerical results and experimentally observ- 



able quantities. In what follows, we begin by briefly 
reviewing selected aspects of the numerical techniques 
used, and how they relate to the physical problem un- 
der study. We then recall the connection between struc- 
ture factors and averaged correlations in random systems, 
and discuss the extraction of effective correlation lengths 
from our numerical data for correlation-function statis- 
tics. Next we exploit the connection between field- and 
correlation function distributions at high temperatures, 
in an attempt to derive approximate analytical expres- 
sions for the latter; such formulae are used in turn, in 
order to compute the corresponding structure factor. A 
short section is dedicated to a reanalysis of the asymp- 
totic behaviour of the widths of correlation-function dis- 
tributions, first presented in Ref. and now comple- 
mented by additional data. In the next section, an ad- 
ditional uniform applied field is considered: free energies 
and uniform magnetizations are calculated on strips of 
both pure and RFIM systems. These quantities are used 
to calculate the corresponding Gibbs free energy which, 
in turn, gives the distribution function for the uniform 
magnetization in zero uniform field. Numerical data are 
then analyzed via finite-size scaling and crossover argu- 
ments. A final section summarizes our work. 



II. NUMERICAL METHODS AND D = 2 RFIM 

We consider strips of a square lattice of ferromagnetic 
Ising spins with nearest-neighbour interaction J = 1, of 
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width 3 < L < 18 sites with periodic boundary condi- 
tions across. The random-field values hi are drawn for 
each site i from the binary distribution: 
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p(hi) = -[6(hi - H ) + S(hi + H a )] 



(1) 



TM methods are used, on long strips of typical length 
L x = 10 6 columns, as described at length in Ref. ||] 
and references therein, to generate representative sam- 
ples of the quenched random fields. Along the strip, we 
calculate correlation functions (as explained in the next 
paragraph), as well as free energies and magnetizations 
(details in Section (v|). 

Here we calculate the disconnected spin-spin correla- 
tion function G(R) = ((TqCTr), between spins on the same 
row (say, row 1), and R columns apart. Related quan- 
tities, such as correlation lengths, are defined with con- 
nected correlations, {cfq<jr) — (c) 2 , in mind; however, for 
the quasi-one-dimensional Ising systems under consider- 
ation (either pure or random) one is always in the param- 
agnetic phase, so the distinction between connected and 
disconnected correlations is unimportant. In Ref. || we 
explained why the ranges of spin-spin distance, temper- 
ature and random-field intensity of most interest for in- 
vestigation by TM methods are, respectively, R/L ~ 1, 
< T < T c o = 2.269 • • • [we take k B = 1], H < 0.5. 
Here we restrict ourselves to high T > 2.0, and rather 
low fields, H < 0.1 — 0.15 . We use a linear binning for 
the histograms of occurrence of G(R); usually the whole 
[— 1, 1] interval of variation of G(R) is divided into 10 3 
bins. 

Since we shall be dealing with probability distribu- 
tions, a word is in order about multifractality. Though 
multifractal behaviour has been found at the critical 
point of random-bond Potts systems [|l0|jll| , the avail- 
able evidence strongly suggests that, off bulk criticality, 
correlation functions behave normally flQl . Thus, in the 
present case we expect that analysis of different moments 
of the probability distribution of G will yield essentially 
the same results. 



III. CORRELATION DECAY 

The properties of correlation functions are usually in- 
corporated into associated correlation lengths, whose ba- 
sic definition is as (minus) the inverse slope of semi- 
logarithmic plots of correlation functions against dis- 
tance. In this view, one assumes both that exponential 
decay can be well-defined at essentially all distances, and 
that a single length is enough to characterize such be- 
haviour. In cases as the present, quenched randomness 
implies that configurational averages must be taken, and 
one must be careful in deciding what quantities are to be 
thus promediated. Recall that, e.g. in neutron scattering 
experiments, the intensity of the magnetic critical scat- 
tering is proportional to the average (over the crystal) of 



the scattering function S(q), which is the Fourier trans- 
form of the correlation function for wave- vector transfer 
q With G R = G(R), and wave vector q in the 

row direction, S becomes 



[S(q)] 



[J dRe lqR G B ] = J dRe 1 * 11 (G R ) 



(2) 



where [• • •] stands for configurational average, and 

(G R ) = J cIGrP(Gr)Gr , (3) 

where P(Gr) is the probability distribution for Gr. The 
last equality in Eq (g) depends only on the assumption 
that P(Gr) is position-independent along the crystal. 

The simplest assumption for P(Gr) that incorporates 
both disorder and exponential decay given by a single 
length £ for all distances is a Gaussian distribution: 



P{Gi 
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/2wA(R) 



-y 2 /2A 2 (R.) 



y=Gj 



-R/e 



(4) 



where distance-dependent widths A(i?) allow for, e.g., 
(disorder-induced) larger uncertainties for larger spin- 
spin separations. However, using Eq. (g) in Eqs. (@)-(|) 
one obtains a width-independent Lorentzian form for the 
average structure factor: 



[S(q)} 



a 2 4- — 



(5) 



This coincides with the standard mean-field result for 
the disordered phase, and is deemed unsatisfactory upon 
comparison with experimental data |7],[ll| . 
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FIG. 1. Normalized histogram P(G) of occurrence of G. 
Strip length L x = 10 6 columns, binwidth 2 x 10~ 3 . Verti- 
cal bars located respectively at: Go (full line), (G) (dashed). 
Shaded region on horizontal axis from (G) — W to (G) + W. 
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We now exhibit our numerical results, and compare 
their implications to those of Eqs. @) and (||). For high 
temperatures and low random-field intensities, as speci- 
fied above, we recall (see also Ref. Q) the following main 
features found for the probability distribution P(G): (i) 
a clearly-identifiable, cusp-like peak, at some G m below 
the zero- field value Go = G(Hq = 0); (ii) a short tail 
below the peak and a long one above it, such that (iii) all 
moments of order > of the distribution are above Go. 
In Figure [lj where the first moment (G) is shown, one 
has G m = 0.278, G = 0.2853, (G) = 0.321; the RMS 
width W = ((G - (G)) 2 ) 1 / 2 = 0.071. 

Therefore, the features depicted in Fig. [I], especially 
the asymmetric cusp, are at variance with the form 
Eq. (^) . We now investigate what effects are carried over 
to the associated correlation lengths. We do so by mim- 
icking the procedure outlined in Eqs. and (^) above: 
first we average over randomness for a given spin-spin 
separation, and then study the variation of the averaged 
quantities over distance. The results are in Fig. ||, where 
our numerical data for (G) are plotted against varying 
R. H — data are also shown for comparison. 
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FIG. 2. Correlation decay along strips of widths L — 5 
and 9. Full lines: H = 0. Points: (G) for H = 0.1. Dashed 
lines: unweighted least-squares fits of Ho = 0.1 data. Vertical 
bars give RMS widths W of corresponding distributions. 

One sees from the respective slopes that, taking into 
account data for R > L, the correlation length for Hq ^ 
would seem to be systematically larger than in zero 
field. This reflects the domain structure into which the 
system breaks down: at short distances the conditional 
probability for a spin to belong to the same domain as 
the one at the origin is larger than for Hq = 0. 

For longer distances, correlation functions start to 



show severe disorder-induced fluctuations, related to the 
crossing of domain walls. Contrary to the zero-field case, 
where temperature-induced domain walls are present but 
the respective sign changes in correlation functions av- 
erage out to give an exponential fall off, here the do- 
main wall configurations are essentially determined by 
the (quenched) accumulated random-field fluctuations. 
At large R, such fluctuations play a very sensitive role 
even for very low random field intensities . One antic- 
ipates problems with defining correlation lengths from 
the corresponding data. The vertical bars in Fig. || show 
that, for fixed strip width L, the width W of the distri- 
bution indeed grows apparently unbounded for increas- 
ing R; though this is related to the crossing of domain 
walls just mentioned, it is also, and predominantly, an in- 
trinsic feature of the quasi-one-dimensional systems used 
here. Thus, inferring two-dimensional behaviour from 
such trends may be risky. However, we now argue that 
in d = 2 one does actually run into problems for large R, 
exactly as inferred above; only, the underlying reasoning 
is subtler. 

In fact, see Ref. and Section [V] below, a different 
analysis of correlation functions, at fixed R/L, strongly 
suggests that the relative widths W = W/(G) grow as 
R, L — > oo in d = 2, approaching a finite limiting value 
G Hq , G ~ 2, k ~ 0.5. This means that, when one con- 
siders the dispersion of ln(G), the signal-to- noise ratio 
becomes of order one for large R, L, and it is this latter 
fact that, in d ~ 2, must compromise attempts to extract 
correlation lengths in such range. 

The effect of the above on fits of neutron-scattering 
data to lineshapes is that, since the latter rely on the 
idea that correlation lengths are always reliable quan- 
tities, they may be off the actual picture in the small- 
wavevector region. 

We now attempt to derive approximate analytical ex- 
pressions for P(Gfl); our ultimate goal is to predict a 
form for [S(q)} from Eqs. (§)- (|). 

IV. DISTRIBUTION OF G FROM FIELD 
DISTRIBUTION (AT HIGH T) 

In this section we use simple scaling ideas to establish 
a quantitative connection between the underlying distri- 
bution of accumulated fields and that of the correlation 
functions themselves. 

We begin by considering a one-dimensional system 
with sites denoted by i = 0,1,..., uniform nearest- 
neighbour interactions K and site-dependent random 
fields hi (both in units of T). For any given specific 
realisation of the fields an exact decimation scheme with 
length scaling factor 6 = 2 can be applied, eliminating 
all odd- numbered sites. The renormalised fields at, and 
coupling between, e. g. spins and 2, are given by: 

K - h' 2 = ho ~ h 2 
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h' + h' 2 = /i 



4if ' = In 



In 



cosh 2 (X + hi) 
cosh 2 (if - fti) 



cosh 2 (if + /ii) cosh 2 (if - /ii) 



cosh 



(6) 



Iterating this procedure n times, one obtains a single 
renormalised bond K connecting sites and R (R = 2"), 
at which the respective rescaled helds are ho, hp. The 
correlation function G(R) = (o _ o <t -r) is therefore: 



G{R) 



a 2K 



cosh(/io + hp) — cosh(/io — hp) 
e 2K cosh(ft,o + hp) + cosh(/io — hp) 



(J) 



For low fields Ho <C 1, one uses cosh(/io — hp)/ cosh(/io • 
/lit) ~ exp(-2h hp) to get 



G(i?) ~ tanh(if + ft, /i R ) 



(8) 



Then, provided also iio *C if, the distribution of G(i?) 
is given by that of X = h hp, since (to lowest order in 
i?o) if' = ilncosh2if is field- independent. One has: 

P(X) = Jdh P(h ) J dhpP(hp)5(X-h hp) . (9) 

At low ii , the scaling equations (||) give h' ~ ho + 
hi tanh2if. Repeated applications of this transforma- 
tion give ho ~ SiLi ^ ^ where £ ~ ln(tanh if ) _1 
is the correlation length at low iio- Then ho (and simi- 
larly hp) is the sum of Af independent variables (Af = R), 
so the individual distributions of ho, hp become (at large 
R, £) Gaussians of width A = HoVAf : 



(10) 



P(h .p) cx exp -(V-R./A) 



For R > £, the same form applies, but because the field 
accumulation under scaling is cut off by the decreasing 
tanh if, the relation for A involves Af ~ £. So, in general, 
Af ~ min (R,(). 

Making ho = s cos 6, hp = s sin9, 



/■00 /-27T 

P(X)(x / dss d0e- s2 / A2 6(X 
Jo Jo 

with the final result 



■sin 2(9) (11) 



P(X) = -£- 2 e-* In(l + i 



21X1 



(12) 



where a is an overall normalization constant. Strictly 
speaking, Eq. ( |l2| ) is the asymptotic reduction of Eq. (11), 
valid for the regimes y <JC 1 (the relevant one for our 
purposes, as shown below) and i/> 1. 
Transforming back to P(G), one sees that the value G m 
for which P(G) is maximum must correspond to X = 0, 
which maximizes P(X). Thus, from Eq. m), 



2/ = -f-|tanh- 1 G ro -tanh- 1 G| 



(13) 



For G close to G„ 



P(G) 



exp 



linearization gives 
-G\ 



-3? I G 



2A 2 [1 - G 2 ] 



In 



|G m - G| 



(14) 



with A 2 = ±A 2 (1 - G^). 

The main feature exhibited by this form is a locally 
symmetric cusp, with infinite slope on cither side, at 
G = G TO . This is expected to carry over to more gen- 
eral contexts, provided that Ho -C 1. Indeed we have 
checked that a similar description (applying approximate 
Migdal-Kadanoff scaling calculations), with the predic- 
tion of a cusp, also applies on strips and in two dimen- 
sions (see Eqs. Jig), ( |l6| ) below, and related discussion). 
A quantitative test of Eq. ( |l4| ) is shown in Figure^, where 
only data for G < G m arc displayed (we shall deal with 
G > G m immediately afterwards). The conditions are 
such that R < £, G 2 m -C 1 (see Fig. ||), so the above 
low-field theory gives A ~ Ho\fR. 
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FIG. 3. Semilogarithmic plots of distribution functions be- 
low Gm, the value at which P(G) peaks (see text). L — 5, 
T = 2.0. Straight line is guide to the eye. 



One sees that exponential decay against \G—G m \/RH 2 
is indeed the dominant behaviour, provided that Ho ^ 
0.15; already for ii = 0.25, small departures show, 
which become more prominent for Ho = 0.5. 

As regards cusp asymmetry, not predicted by Eq. (|l4|) , 
we have found that although data for G > G m still fall ex- 
ponentially for small Ho, they do not collapse when plot- 
ted against |G — G m \/ RHq. This is because the mutual 
reinforcement, between ferromagnetic spin-spin interac- 
tions and accumulated field fluctuations (responsible for 
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the long forward tail Q), is left out by the approximation 
above Eq. (H), namely that K' is £fo-independent. 

Before calculating [S{q)\ from Eq. (|lj), we recall that 
Eqs. (d)-(H) are normally required for bulk systems, thus 
one must work out an approximate scheme to go from 
the d = 1 regime of Eqs. (§)-(|l|) to d = 2. We have 
done so via a Midgal-Kadanoff rescaling transformation 
at T ~ T c q. As a consequence of the similarity of the 
corresponding recursion relations, to those for one dimen- 
sion, one ends up, after m scalings such that 2 m = R, 
with a result very similar to Eq. (fsl) : 



G(R) = tanh(i^ + h h R ) 



(15) 



where again one assumes low fields, ho, For large 
R these have Gaussian distributions of width A/j deter- 
mined by the eigenvalue A of the low-field scaling trans- 
formation of Hq. For T near T c o, A# cx i? M Hq where 
fi = In A/ In 6 (b = 2). Further, 



K~K r 



R 

1 



(16) 



Since Eqs. j9|)-(|l2]) still apply, provided A is replaced by 
A#, one gets the dominant contribution to the scattering 
function as: 



[S{q)} cx Re 



dRe iqR e 



R/S 



1 



CAl) 



(17) 



where C is a constant of order unity. This can be trans- 
formed into: 



[S(q)} cx 



■CH^(2fi)\ Re 



i2 _ 2 , ^„ ov . w ._^-.g 



1 



(18) 



If one assumes the form A = Hq\/R, given for R < £ 
in one dimension (see above and below Eq. (|Io|)), and 
also by the Migdal-Kadanoff scheme in d = 2, then 
/i = 1/2 and Eq. (H) predicts the lineshape to be 
Lorentzian plus Lorentzian-squared, the mean-field form 
found when the disconnected contribution is taken into 
account ||M|l2|-|l3] . On the other hand, if ones goes by 
the saturation behaviour predicted in Ref. and Sec- 
tion [y| below, and by the scaling approaches if R > £, 
then the result is C = 0, corresponding to a single 
Lorentzian in Eq. (|lq). 

Though cither of these final predictions is certainly 
open to challenge, in view of the number and severity 
of approximations involved in the course of their deriva- 
tion, it is expected that the procedure described above 
will serve as a rough guide to attempts at connecting 
basic microscopic features (such as fluctuations of accu- 
mulated fields) to observable quantities, e.g. scattering 
functions. 



V. WIDTHS OF G-DISTRIBUTION 

In Ref. |^| we studied the variation of the RMS relative 
width W of the probability distribution of correlation 
functions, against field intensity and strip width, for fixed 
R/L, high temperatures and small H . We proposed the 
scaling form 



W = Hq /(LHq) 



(19) 



J c 0, 



good data 
can indeed 



and showed that, for R/L = 1, T 
collapse of y = In [W Hq K ] against x 
be obtained with k ~ 0.43 — 0.50 and u ~ 0.8. We used 
L < 15 and scanned < x ^ 2.8; keeping k = 0.45 
and u = 0.8, we found for x > 1 a satisfactory fit given 
by y = —0.3 — 5.3 exp(— 1.57a;), which would imply an 
exponential saturation of the scaled width W hy K as x — > 
oo, with a limiting value exp(— 0.3) = 0.83. 

In Figure |[ we display again the data of Ref. || , plus 
additional data for L = 15 and 18, which enabled us to 
explore larger values of x (x < 4.0) while still keeping 
to relatively low Hq. We then reanalyzed our full set of 
data, with the results that (i) we managed an excellent 
fit to the whole interval < x < 4 by a single expression 



u : 2.2 { in [ — 



0.6 



0.4025 



(20) 



(in which y and x involve the same values of the expo- 
nents k, u as earlier) and that (ii) this new fit, while 
still predicting saturation for x ^> 1, implies that in 
the approach to two dimensions x ^> 1, convergence of 
W hq* 1 is power-law-like, giving a limiting scaled width 
W tiQ K = 0.6~ 2 - 2 exp(-0.4025) = 2.06. 
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4. Semi-logarithmic scaling plot of RMS relative 
W Hq against LHq. Curve is fitting spline, given 
l20|). 
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VI. MAGNETIZATIONS 

In this section we examine the scaling properties of the 
uniform magnetization on strips of the d = 2 RFIM. For 
convenience we shall always keep T = T c q. 

We first outline our method, which involves a gener- 
alised Legendre transformation. Consider the Hamilto- 
nian 



(21) 



where Ui are Ising spins, and Tig includes all interactions 
except that of the spins with the uniform field h. One 
has, for the corresponding partition function Z{h): 



Z{h) 
Z(0) 



(3AIh 



A I 



(22) 



where Pq(M) is the probability of occurrence of the value 
M for the magnetization, in zero uniform field. Assum- 
ing a system with N 1 spins, with f{h) = negative 
free energy per site in units of T, and m = (3AI/N, 



o N(f(h)-f(0)) 



N / dm Pq (m) e 



N mh 



(23) 



In order for extensivity to be satisfied, one must have 
7VPo(m) = exp N g(m), where g{m) is intensive and de- 
termined by 



N(f(h)-f(a)) _ / dme N{9{m)+rnh) 



(24) 



Assuming the usual sharp-peaked distribution around a 
thermodynamically averaged value m, one sees that 



f(h)- f(0) = g(m) + mh + O 



In A 



(25) 



with (dg/dm)jn = —h. That is, g is the standard Gibbs 
free energy per site. Substituting back in Eq. (^3|), one 
gets: 



Pa(M) = expA<?(m) , 



(26) 



where terms of C(ln N/N) have again been neglected. 
Eq. (26), with g(m) given through Eq. (p5|), is the nat- 



ural starting point to study magnetization distributions 
by TM methods. Indeed, though one can get the thermo- 
dynamically averaged exact values of all moments of the 
distribution via TM (ll],[l(|, the distribution itself is not 
given directly. This contrasts with Monte-Carlo meth- 
ods, which inherently incorporate readily-observable fluc- 
tuations around equilibrium, and have been widely used 
to study magnetization distributions at criticality, both 
in hypercubic geometries and on planar lattices with 
various aspect ratios |lq] . 

Recall that, on strips of width L and length L x , 
L x ^> L such as is the case here, the aspect ratio is 



essentially infinite, therefore Pq(M) will be Gaussian, at 
least for pure systems fl6]jL8l| . Our purpose (as shown 
below) is to compare pure- and RFIM- results and ex- 
plain their mutual differences, by using general theory 
of RF systems [|2|.||.[L9] ^l) coupled with finite-size scaling 
(FSS) || . 



A. Pure Ising systems 

We start by illustrating the properties of g{m) for pure 
Ising spins. One calculates f(h), f(0), m(h) in Eq. ( gg ) 
by standard numerical methods |23|] : the first two by 
isolating the largest eigenvalue Aq of the TM and us- 
ing / = L _1 InAo (which is tantamount to assuming 
L x — > oo; more on this below), the third by calculating 
derivatives of / relative to h. The latter is done here by 
perturbation theory |l5|,[l6|,^4],^5| , both for better numer- 
ical accuracy and because an adapted procedure proves 
convenient when dealing with the RF case, where sam- 
ples over disorder must be accumulated. 
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FIG. 5. Scaling plots of magnetization and excess free 
energy for pure Ising systems at criticality. Strip widths L = 4 
(circles), 8 (crosses), 12 (triangles), 16 (squares). Normalized 
magnetizations (= T c df/dh) are used to avoid superposition 
of plots. 



At t = (T - T c )/T c = 0, FSS |22) gives for the excess 
free energy: Af(h, L) = f(t = 0, h, L) - /(0, 0, L) _ 
h 1+1 / s FiLh 1 ^), with 8 = 15, y h = 15/8. In Fig. ~ 
we show scaling plots of Af(h) /i _ ( 1 + 1 /' 5 ) and mhr 1 
against Lh l ^ Vh . For low fields (h < L~ Vh ), the slopes of 
both logarithmic plots are given by the (finite-size) initial 
susceptibility exponent 7/z/ — 7/4, as a consequence of 
the scaling relation y^(l — 1/5) = j/v. 
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For t non-zero, but still for low fields, one generally 
expects Af(h) — a(t,L)W l whence m — a(t, L) /i h^ 1 , 
g = Af — mh = a(t,L)(l — /i)/i M , implying g ~ 
afaLyV^-Vm^^-V. Subcases are: 

(i) t = 0, L = oo : Af ~ so 5 - m 1+s ; 

(ii) i small, L = oo, 1 > t > to 1 /* 3 : A/ = a(t) h 2 , a(t) = 
5X(t)~*~ 7 , so S ~iV; 



(iii) t = 0, L finite, 1 < L < 
o(L) ft 2 , a(L) = ~ £ 7/v , so ff ~ L'^ ; 



I M- 



A/ = 




iQ- 4 icr 3 i(r 2 o.i 1 

2 

m 

FIG. 6. Negative Gibbs free energy —g = — (A/ — m/i) 
times I/ 7 /" against m 2 for pure Ising systems at criticality. For 
key to symbols, see caption to Fig. H. Straight line has unitary 
slope and is a guide to the eye. Normalized magnetizations 
on horizontal axis only. 

Case (iii) is depicted in Fig.|(]. One sees that — g Lp/ V ~ 
m 2 as far as m ~ 0.6, which (through Eq. (p6|)) is con- 
sistent with the Gaussian behaviour predicted for Pq(M) 
in this case. Close to m = 1 scaling breaks down, and 
the deviation from saturation magnetization must follow 
a single-spin-flip picture, e = 1 — m ~ exp(— 2/T c ). The 
effects of this on g can be worked out from a high-field 
expansion, in which TLq of Eq. (|2l|) is taken as a pertur- 
bation on the field term h^ li tjj |26[. The result is: 



ds 



so 



W- 

2 £ 



(27) 



where s = \ (ln2-8/T c ) = -1.41617.. . . 

Figure shows that Eq. (|27| ) is in excellent agreement 
with numerics already for (1/2) ln(l/e) ~ 2.4 (e ~ 10 -2 ). 
This provides a rigorous check of our analytic and numer- 
ical procedures. 
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(l/2)ln(l/e) 
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FIG. 7. Derivative of Gibbs free energy relative to magne- 
tization deviation e against (1/2) ln(l/e). Points: numerically 
calculated derivatives from L — 16 data for /, m (L — 8 al- 
ready gives results indistinguishable from those displayed). 
Straight line: first two terms on RHS of Eq. (Ej). 

Returning to the connection between g and magneti- 
zation distributions, we first note that, although we are 
using / = L _1 InAo which holds only for strip length 
L x — > oo, the number N = L X L of spins in Eq. ( pt^ ) 
implies a finite, though possibly very long, strip for the 
Equation to be of practical use. An estimate of the er- 
ror implied by using infinite-strip free energies instead 
of their fully-finite system counterparts can be obtained 
by referring to Table 2 of Ref. |16|, where it is shown 
that for systems with aspect ratio a = L x /L = 100, 
the corresponding value of the Binder parameter Q = 
(M 2 ) 2 /(M 4 ) is ~ 2% off its a -> oo limit, 1/3. As typi- 
cal widths used here and, especially in the next section, 
are in the range L < 20, and assuming that errors in 
the distribution and in its calculated moments are of the 
same order, it follows that using the infinite-strip expres- 
sion for / implies deviations in Pq(M) smaller than 2% 
for L x > 2, 000. The advantage of this procedure is that 
L x can be seen essentially as a free parameter, i.e. not 
connected to an actual number of iterations along the 
strip. For our purposes here, we shall always be com- 
paring pure-system results (where L x is fictitious in the 
sense just described) with those obtained on RF strips of 
the same width at the same temperature, where sampling 
over disorder typically necessitates an actual L x ~ 10 5 ; 
thus, equating the values of L x in both systems is both 
correct as far as comparisons are concerned and, given the 
lengths required for adequate sampling over randomness, 
fully within acceptable error margins for the description 
of pure systems. 



From Eq. (p6|), for case (iii) where —gL 1 l v 



9o m 



7 



<7o — 10 1 as shown in Fig. 



P {M) = exp 



-go 



(28) 



where 2 — r\ = 7/z/ was used. Therefore the width of the 
Gaussian distribution is W ~ g^ 1 '" 2 L^- 1- ^/ 2 / yfi^. 

B. RFIM 



We now include the term ^ /ijCj in 7Yo of Eq. (21), 
with the local fields /i, distributed according to Eq. (|l]). 

We first consider the application of FSS to the RFIM 
in zero uniform field. For bulk systems, theory pre- 
dicts that the scaling behaviour of the RFIM depends on 
Ho|£| where Hq is the random-field intensity, t = (T — 
T C (H ))/T C {H Q ) is a reduced temperature §,|,(jj|-fl|, 
and Q the crossover exponent is <j) = 7i the pure Ising 
susceptibility exponent. For d > 2, T c (Hq) is the field- 
dependent temperature at which a sharp transition still 
occurs; in d = 2 the dominant terms still depend on the 
same combination, where now ||l[] "T c (iJo)" denotes a 
pseudo-critical temperature marking, e.g., the location 
of the rounded specific-heat peak. In d = 2, specific 
heat (2(]] and neutron-scattering |^lj data are in good 
agreement both with the choice of scaling variable as 
above, and with the exactly known 7 = 7/4. For the 
excess free energy A'/ = f(t, Ho) — f(t, 0) in two dimen- 
sions, an additive logarithmic correction arises 



A'/ = A*t 2 In H + H^^^it H Q 



(29) 



In Ref . H , we showed that the appropriate FSS variable 
for the description of correlation functions in finite RF 



LH 



2v/4> 



While the second term 



systems at t — is x ; 

on the RHS of Eq. (^9j) is in that way taken care of, the 
logarithm needs separate consideration. On the basis of 
renormalization-group arguments, in which L~ l is seen 
as an additional relevant field p8) , one realizes that the 
steps leading to the appearance of the t 2 In Hq term in 
Eq. ( |29| ) also apply here. Indeed, since the respective 
eigenvalue [yx = 1 in that case, y^-i = 1 here] divides the 
dimensionality d = 2 po| , ^7| ), a corresponding scenario 
obtains at t = and L ^> 0, when L~ x l v is substituted 
for t. Therefore, we assume: 



A'f(t = 0,L,H o ) = AL- 2 lnH 



Hl du/ ^{LH { 



) 



(30) 



In Fig. pi where T = T c (thus a small, ^"dependent 
shift in U T C (H )" has been neglected, which 

should not matter much for low RF intensities), are dis- 
played results of a numerical test of Eq. ( |30| ) , both with 
and without the logarithmic term. 

We have found fits of a quality similar to that shown 
in Fig. H (b), where A — 10~ 4 , for a wide range 10~ 5 <■ 



A < 10 2 along which the \ 2 estimator remains approxi- 
mately constant. At large Hq L^l 2v , however, the fits de- 
teriorate noticeably (not obvious from Fig. || (b) , because 
of the large vertical scale) , no doubt owing to the incip- 
ient breakdown of the small-RF regime (where, e.g. the 
.ffo-dependent shift in "T c (i? )" is no longer negligible). 
Comparison with experimental data e.g. from Ref. po| is 
not straightforward, as transforming from bulk scaling, 
Eq. (|2^), to FSS, Eq. (|30|), may involve numerical factors 
not immediately available. 



< 




a 



10 '■ 10 : - 0-1 
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FIG. 8. A/ s = (A'/(L,// )- J 4 J L- 2 lni?o)-ffo~ 2d " / ' /> P lotted 
against H L 4 ' /2,y 1 4> = 7/4, v = 1 . (a): A = 0. Bottom to 
top: L = 4, 6, 8, 10, and 12. (b): A = 10 -4 , same notation. 



Moving on towards incorporating both RF and uniform 
field effects, we again neglect the -ffo-dependent shift in 
u T c (Ho)" and make T — T c o- In the presence of sev- 
eral relevant fields u\,U2 ■ ■ ■ with respective scaling pow- 
ers yi,y2 1 ■ ■ ., the singular part of the free energy scales 

f(ui,u 2 , ■■■) = M dM F ( , U2 , , , U \ , ■ 

(31) 

Using u\ = Hq, U2 = L , U3 = h, one has in the case 
(d = 2): y 1 = 2^/0 - 8/7, y 2 = 1, y 3 = y h = 15/8, 
therefore 

f(H Q , L, h) = Hl dv/,p F(L Ho V/ *, h H- 2vVh/ *) . (32) 

Possible In Hq corrections in the manner of Eq. ([50|) have 
been omitted, since our interest will focus on the calcu- 
lation of the Gibbs free energy, which in the case de- 
pends on A"/ = f(H ,L,h) - f(H ,L,h = 0) (see 
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Eqs. (|21|)-(|25|)); we are thus assuming that, at least for 
small enough h, the logarithmic terms cancel in the sub- 
traction. 

Similarly to the pure case but always at t = and 
L- 1 - 0, we investigate the small- h regime, in which 
one expects A" f = a{L,H Q )hP. From Eq. @, this 
implies 



A"/ = hP Hq 



2v{d- 



(33) 



By assuming, as Hq —>■ 0, a power-law dependence 
Fi(x) ~ x , and demanding that, in this limit, (i) the 
ifo - dependence of A"/ must vanish and (ii) the form 
h 2 L J / U be reobtained, one gets fi = 2, t — 7/4. There- 
fore, one has generally for small h < L~ Vh : 



A"/ 



h 



(34) 
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FIG. 9. Negative Gibbs free energies: (a) raw data and (b) 
scaled, against squared uniform magnetization. Linear scales 
used on axes, in order to underline spread of raw data. In (a) , 
data for same (L, Ho) are joined by dashed lines. Strip widths: 
L = 4 (circles), 6 (crosses), 8 (triangles) and 10 (squares). 
RF intensities: Ho = 0.3, 0.4, 0.5 (L = 4); for L = 6, 8, 10, 
Ho = 0.2, 0.3, 0.4. For fixed L, Ho increases from bottom to 
top curves. In (b), plots of —LHo~ 6 g collapse well onto a 
straight line against m 2 (Eq. (^6|) with z = 1) up to m 2 ~ 0.4. 
Straight line in (b) is guide to the eye only, and has slope 
0.275 . 



Considering now a non-vanishing, but still small, RF, 
one may assume crossover to a new power-law form for 
F x {x), Fx{x) ~ x z . In this regime [LHf l4> > 1), 



A"/ = L z H, 



2vz/<p-2 ,2 
'* 



whe nce a (L, Hp) = L z Hq VZ ^ 2 , yielding (see Subsec- 



tion 



VIA) 



g(L, H , m) = - [a(L, H )] 1 m 2 



(36) 



Our data for g(L, Hq, m), displayed in Fig. g, are consis- 
tent with z = 1, that is, g(L, H , m) 
One then has, using Eq. (|2q), 



6 / 7 , 



P (M)=exp - gi L x H^ /7 m 2 



(37) 



with gi ~ 0.28 from the slope of the straight line in Fig. 0. 
Therefore the distribution is still Gaussian, with a width 
W ~ g x X ^ 2 H Q ' I \fL~z. Comparison with a correspond- 
ing pure system, see Eq. ( p8| ) and the arguments in the 
paragraph preceding it, gives: 



W{L,L X ,H Q ) 
W(L,L x ,0) 



(lhT) 



-3/8 



(38) 



showing again that the FSS variable x = L Hq V ^ is the 
relevant one. For x 3> 1 where RFIM behaviour sets in, 
one sees that distribution widths are smaller for RFIM 
than in zero field. 



(35) 



VII. CONCLUSIONS 

We have used TM methods to calculate spin-spin cor- 
relation functions, Helmholtz free energies and magneti- 
zations on long strips of width L = 3—18 sites of the 
two-dimensional RFIM, close to the zero-field bulk criti- 
cal temperature. 

Through analysis of the probability distributions of 
correlation functions for varying spin-spin distances R, 
we have shown that fits to exponential decay of averaged 
values against R (for R not too large) give rise to effec- 
tive correlation lengths larger than in zero field. This is 
because of the reinforcement of correlations within do- 
mains. At longer distances (i.e. across many domain 
walls, R/L ^S> 1), fits of exponential decay become unre- 
liable, thus compromising definitions of effective correla- 
tion lengths. 

We have worked out explicit connections between field- 
and correlation function distributions at high tempera- 
tures, yielding approximate analytical expressions for the 
latter. Such expressions account well for trends found in 
numerical data, namely the existence of peaked cusps and 
the functional dependence, on R and field intensity Hq, of 
data below the peak; above the peak, although agreement 
with numerics is not good, we have pinpointed that the 
responsibilty for this lies in a truncation in our approx- 
imate scaling scheme, which decouples scaled nearest- 
neighbour interactions from the random field. We have 
discussed the use of analytical expressions, such as the 
ones found here, for computation of the corresponding 







structure factor. Though results as they stand are far 
from conclusive, we have established a rough guide to 
attempts at connecting basic microscopic features, such 
as fluctuations of accumulated fields, to experimentally 
observable quantities, e.g. scattering functions. 

We have reanalyzed the asymptotic behaviour of the 
relative widths of correlation-function distributions, first 
presented in Ref. and now complemented by addi- 
tional data. While our earlier analysis seemed to point 
towards exponential saturation, the new set of data 
shows that, for fixed R/L = 1, the fractional widths 
of correlation-function distributions behave consistently 
with asymptotic power-law saturation, i.e. depending on 
L~ 22 , see Eq. (p0|). The scaling variabls remain as given 
previously. 

Considering a uniform applied field h, we have derived 
a connection between Hclmholtz free energy f[h), uni- 
form magnetization m(h), the Gibbs free energy g(rn), 
and the distribution function for the uniform magnetiza- 
tion in zero uniform field, Pq(to), which is in principle 
applicable to any finite system. By working at the bulk 
zero-field critical temperature T c , we have illustrated 
our approach by showing that, for strips, one indeed gets 
a Gaussian distribution |l6|,[l8| for m not very close to 
saturation. Near m = 1, where scaling breaks down and 
a single-spin- flip picture holds, a perturbation expansion 
accounts for the properties of g(m). Still at T c o, now 
in non-zero random field, we have found from finite-size 
scaling and crossover arguments, coupled with numerical 
data, that for strip geometries, Po(m) is still Gaussian, 
and its width varies against (non- vanishing, but small) 

— 3/7 

random-field intensity Hq as H . This is again valid 
far from saturation (typically, for m 2 < 0.4, see Fig. ||). 
The ratio between the width of Po(m) and the width of 
the corresponding distribution for a strip of same length 
and width in zero field varies as (L H^ 7 )~ 3 ^ 8 . 

We expect that at least some of the features discussed 
here, for distributions of correlation functions and mag- 
netizations on strips, translate also for other geome- 
tries. Considering, for instance, square systems: docs 
the non-trivial form of Po(m) at bulk criticality in zero 
field [[l6|-|Is|] evolve into a corresponding shape for H a ^ 
which depends on the variable (L H^ 7 )~ 3 / 8 as here? 

Finally, recalling possible connections with experi- 
ment: given that the description of correlation decay via 
effective correlation lengths runs into difficulties at long 
distances, perhaps this compromises naive fits of neutron- 
scattering data in the small-wavevector region. 
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